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ABSTRACT 

We analyze the first cosmological simulations that recover the fragmentation of high-redshift 
galactic discs driven by cold streams. The fragmentation is recovered owing to an AMR 
resolution better than 70 pc with cooling below lO^K. We study three typical star-forming 
galaxies in haloes of ^ 5 x IO^^Mq at z ~ 2.3, when they were not undergoing a major 
merger The steady gas supply by cold streams leads to gravitationally unstable, turbulent 
discs, which fragment into giant clumps and transient features on a dynamical timescale. The 
disc clumps are not associated with dark-matter haloes. The clumpy discs are self-regulated 
by gravity in a marginaly unstable state. Clump migration and angular-momentum transfer 
on an orbital timescale help the growth of a central bulge with a mass comparable to the 
disc. The continuous gas input keeps the system of clumpy disc and bulge in a near steady 
state for several Gyr The average star-formation rate, much of which occurs in the clumps, 
follows the gas accretion rate of ^ Ah M^yr'^ . The simulated galaxies resemble in many 
ways the observed star-forming galaxies at high redshift. Their properties are consistent 
with the simple theoretical framework presented in lDekel. Sari & Ceverinol (l2009l DSC). In 
particular, a two-component analysis reveals that the simulated discs are indeed marginally 
unstable, and the time evolution confirms the robustness of the clumpy configuration in a 
cosmological steady state. By z ^ 1, the simulated systems are stabilized by a dominant 
stellar spheroid, de monstrating the process of "morphological quenching" of star formation 
dMartig et alj|2009l) . We demonstrate that the disc fragmentation is not a numerical artifact 
once the Jeans length is kept larger than ^ 7 resolution elements, i.e. beyond the standard 
Truelove criterion. 

Key words: cosmology — galaxies: evolution — galaxies; formation — galaxies: kinematics 
and dynamics — galaxies: spiral — stars: formation 
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5^ , 1 INTRODUCTION 



giant clumps ( Cowie et al. 19951: Ivan den Berghetaill 19961: 
lElmegreen et al.1 l2005bi l2007h . Although some of these 
irregular morphologies may be associated with mergers or 
galaxy interactions, the morphology and kinematics of many 
of the SFGs ind icates that they are rotationally supported, 
exten d ed discs (|Genzel et al.l 120061: iForster Schreiber et al.l 
120061: iGenzel et all |2008|) . It seems that most SFGs are 
vjt^iiAi^i ai. | | z.wuo|, l o tmr^ tj^L ai. | | z.uuoh | j^aw i^i ai. | incompatible with being Ongoing ma) Or mergers Or remnants 

I Forster Schreiber et alj|2009 | ) . Both observation and theory ^.r„.rc J^L^iC^.t ^il "FT^ iBournaud et aP 



Observations at high 
1 — 3 have opened 
active period of galaxy 
Adelberger et al. 2004 



redshifts in the range z ~ 
a new wind ow into the most 
formation jSteidel et al.l 119991: 
Daddi et al.l l2004t iGenzel et al.1 



2006; F orster Schreiber et al J 120061: Eknegreen et al.l 



Genzel et all 120081: IStark et al.l 120081: iLaw et al 



2006; 



2009 



indicate that this was a rather violent process. Continuous, 
intense, cold gas streams along the cosmic web fueled 
dense gaseous discs and induced high s tar formation rate s 
(SFR) on the order of - 100 Af© yr"! dPekel et alj|2009l) . 
much higher than the few M© yr^^ in today's galaxies. We 
refer to these galaxies as SFGs, standing for "Star-Forming 
Galaxies" and for "Stream-Fed Galaxies". In these extreme 
conditions, the images of most SFGs are irregular with 
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[maiormei 
of suc h mergers dShapiro et akl 120081; . 
200^ Dekel et al.l 120091; iBtournaud & ElmegreenI 120091; 
iDekel. Sari & Ceverinol l2009l). th ough there are counter 



examples (e.g.. [Robertson & Bullo ck 2008). This indicates 



a paradigm shift in our understanding of galaxy formation. 

The high-redshift galactic discs are very different 
from the low-redshift disc galaxies. In a typical z ~ 
2 disc, 10 — 40% of its UV-restframe light is emitted 
from a few clumps of a characteristic size ~ 1 kpc, 
and with a hi gh SF R of tens of M© yr~ ^ each 



( Elmegreen et alj 2004al ; Elmegreen & ElmegreenI 12005 at 



IForster Schreiber et al.l 120061; IGenzel et al.l l2008l) . These 



clumps represent giant star-forming regions, much larger 
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than the star-forming complexes in local galaxies. This 
pronounced clumpy morphology at high redshift is not 
limited to the UV as the main high-redshift clumps are 
also seen in optical rest-frame emission (iGenzel et al..i2008i' 
iForster Schreiber et a n 120091) . It is not a bandshift effect 
in the sense that the UV images of clumpy low-redshift 
galaxies would not appear so clumpy if put at high redshift 
and o bserved with limited resolution and signal-to-noise 
ratio (lElmegreen et al.l [2009h . The clumpy galaxies where 
first dubbed "clump-clusters" galaxies and "ch ain" galaxies 
when viewed face on or edge on respectively dCowie et al.l 
Il995t lElmegreen et all l2004bl) . Most of the massive ones 
are thick rotating discs, with high velocity dispersions of 
tr = 20 — SOkms^^ (one dimen sional) and rotatio n to 
dispersion ratio of V/a ^ 1 — 7 (ICresci et al l l2009h . as 
opposed to today's thin discs of cr ~ 10 km s~^. The smaller 
galaxies contain a l a rge fr a ction o f dispersion-dominated 
system s dLaw et al.l l2007l 120091: iForster Schreiber et al.1 
|2009|) . Estimates of the total gas fraction in SFGs, based on 
CO measure ments, range from 0.2 to 0.8, with an average o f 
- 0.4 - 0.6 (iTacconi et alJl2008t iDaddi et al.ll2008L l2009h . 
systematically higher than the typical gas fraction of '--^ 0.1 
in today's discs. This gas consists of large molecular cloud 
complexes that host the observed star-forming regions. The 
typical age of the stellar populations in these clumps are 
crud ely estimated to range from one to several hundreds Myr 
( lElme green et al. 2009; Forster Schre iber et al. 2009), on the 
order of ten dynamical times, indicating that the clumps may 
plausibly survive for such durations. 

A substantial 
fractio n of these high-z clumpy discs show a centr al stellar 
bulge (IGenzel et al.l 120081: lElmegreen et al.l l2008a) . Other 
observations indicate that a large population of compact, 
massive spheroids with suppress ed SFR is already present 
at those redshifts dKr iek et al. 120061 : Ivan Dokkum et al] 
l2008HKriek etani20()9l) . This means that a well-developed 
Red Sequence was already in pla ce at z ~ 2. While 



major mergers do lead to spheroids (iRobertson et al. 120061: 
ICox etal.ll2006t iDekel & Coxll2006t iHopkins et al.ll2006h . 



the major-merger rate is not enough for producing spheroids 
in the observed abundanc e ( iDekel et al.l 120091 : iGenel et al] 
120081: IStewart et alj l2009l and references therein), so an 
additional and complementary mechanism must have been 
working efficiently at high redshifts. 
The 

gravitational fragmentation of gas-rich and turbulent galactic 
discs into giant clumps, and their subsequent migration 
into a central bulge, have been addresse d before in the 
simphfied context of isolated galaxie s ("Noguchi' '1999"; 
Elmegreen et al. 2005bl: Bournaud et a 12007. iGenzel et al.i 



20081: lElmegreen et ali lBoSal). Accordin g to the standard 
Toomre instability analysis Tibomrell 19641) , the material in a 
rotating disc becomes unstable to local gravitational collapse 
if its surface density is sufficiently high for its self-gravity 
to overcome the forces induced by differential rotation and 
velocity dispersion that work against the collapse. If the disc 
is maintained in a marginally unstable state with relatively 
high velocity dispersion, the clumps are big. Then, the 
timescale for their migration to the centre by gravitational 
interactions is short, on the order of ten disc dynamical 
times, or a couple of orbital periods at the disc radius. 

The clumps and the other transient perturbations in the 



disc drive further inflow of disc mass toward the centre 
through angular-momentum transport outward. Together 
they lead to a rapid growth of a bulge in about 0.5 Gyr, at 
the expense of the draining and stabilization of the disc. 

The large fraction of clumpy massive galaxie s at z ^ 2 
(e.g., lElmegreen et al.l 120071: iTacconi et al.l l2008l) suggests 
that the clumpy phase should last for a few Gigayears. 
DSC have argued that when a galaxy is part of the 
high-redshift cosmological environment, it can naturally be 
in a long-term steady state typically consisting of a clumpy 
disc and a comparable bulge. Massive galaxies at that 
epoch are fed by continuous, intense col d gas streams that 
follow the filaments of the cosmic web (Keres et al. '20()5t 
Dekel & BirnboimI 120061: lOcvirk et al.i i2008i: .Dekel et akl 
20091) . The fresh supplv of cold gas refills the disc as 
it is being drained into the bulge, and keeps the surface 
density sufficiently high for disc instability and new clump 
formation as long as the bulge is not too massive. On the 
other hand, if a large fraction of the incoming streams 
is already in clumps, they can stabilize the disc against 
fragmentation in two complementary ways: by rapid direct 
buildup of the bulge via mergers and by generating high 
velocity dispersion in the disc. This led DSC to suggest that 
a low level of dumpiness in the streams is a key factor in 
maintaining a clumpy disc in steady state. 

This long-term clumpy phase can be naturally addressed 
by cosmological, hydrodynamical simulations, though the 
numerical challenges are not trivial. The dynamical range 
must explore the cosmological streams on scales of hundreds 
of kpc while properly resolving the gas physics in the 
disc on scales of tens of parsecs. These difficulties 
limited the earlier simulatio ns of disc fragmentation to 
idealized discs in isolation dNoguchil 119991: llmmeh et alj 
2004allH: i Bournaud et al.' '20071; 'Elmegreen et al.' l2008al; 
Bournaud & Elmegreen 2009; Bournaud et al. 2009). They 
successfully 

reproduced the Toomre instability, the subsequent migration 
of clumps and the formation of a classical bulge. However, 
in the absence of fresh cosmological gas supply, these 
phenomena were limited to one episode of instability and 
disc evacuation, and were not capable of capturing the 
long-term cosmological steady state. On the other hand, 
the earlier simulations that did follow the cosmological gas 
supply did not have the necessary resolution with proper 
gas physics to reproduce the gravita ti onal f ragmentation 
of the discs ([Springel & Hernquisl _ 2003 ; jKeres et al 
20051: iGovernatoet al.ll2007l:IOcvirk et al.ll2008l: iDekel et al 



2009) 



A key ingredient for recovering the fragmentation is 
allowing the gas to cool to temperatures below 10^ K 
and reach densities weU above 1 atom per cc, typical 
of gas in giant molecular clouds. If the gas is kept at 
higher temperatures, the pressure prevents local gravitational 
collapse. This may be fine in the modeling of present-day 
disks that do not suffer from wild gravitational i nstability 
dSpringel & Hernquistl2003l: IRobertson et al.ll2004 . but it is 
not adequate when simulating the high-redshift gaseous and 
continuously fed disks that are likely to develop wild disk 
instabilities. 

Even with state-of-the-art resolution of ~ 30 pc, where 
the gas can cool to ~ 100 K and reach n 100 cm~'^ 
dCeverino & Klvpinll2009l; lAgertz et al]|2009bl) . flie thermal 
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Jeans scale in the dense gas is below the resolution 
scale. In other words, there is no pressure to support the 
gas against gravitational collapse on the res olution scale, 
so the gas tends to fragm ent artificially (iTruelove et al.l 
119971: iBate &Burker^ll997l) . In order to prevent artificial 
fragmentation, the common procedure is to introduce a 
pressure floor ( Machacek et al.l200ll : [Robertson & Kravtsovl 



l2008HAgertz et al.ll2009air but it has to be applied in a way 
that would not prevent real fragmentation on the turbulence 
Jeans scale when it is only a few times the resolution scale. 

In this paper we present the results from three different 
high-resolution cosmological simulations, zooming in on 
three nearly random galaxies with haloes of ~ (1 — 2) x 
IO^^Mq at z = 1. They have reached - 5 x lO^Mg 
at z = 2.3, where we perform a detailed analysis. In 
these simulations, we do allow cooling to temperatures 
well below 10*K and thus allow for real physical disc 
fragmentation, while applying the pressure floor necessary 
for preventing artificial fragmentation. We find that they all 
show clumpy discs similar to the observations and according 
to theory. In |2] we summarize the theoretical framework 
laid out in DSC. In |3]we describe the simulation method. 
In f|4] we address the global properties of the simulated 
galaxies and verify that they obey the observed scaling 
relations at high redshift. In f|5]we study in detail the three 
simulated galaxies at z = 2.3 where we focus on the 
disc fragmentation and bulge formation, in comparison with 
theory and observations. i)5.2| brings projected images of 
the different galaxy components, compares simulated 
and observed images, describes the properties of the 
clumpy discs and bulges in comparison to the theory as 
laid out in DSC, S|53] refers to a two-component instability 
analysis, ^5.6| follows the formation of disc clumps and their 
migration to the bulge, and i)5.7| highlights the cosmological 
steady state of clumpy discs, and the possible eventual 
stabilization by a massive spheroid. In Sj6] we test the 
robustness of our results to artificial fragmentation. Finally, 
in 0we summarize our results and discuss them. 



2 THEORETICAL FRAMEWORK 

In a companion paper, DSC have laid out a simple theoretical 
framework for clumpy discs driven by cosmological cold 
streams at high-redshift. We summarize the theory here to 
help us interpret the simulation results as we present them. 
Note that the theoretical expectations did not affect the 
numerical experiments themselves in any way. 

Well before z ^ 1, the high accretion rate is expected 
to keep the discs gas rich, which validates a simplified 
one-component analysis as a first crude approximation. 
According to t he standard Toomre inst ability analysis 
dToomrd Il964t iBinnev & Tremaind l2008l Chapter 6) a 
rotating disc becomes unstable to axi-symmetric modes 
once the local gravity overcomes both differential rotation 
and pressure due to turbulence or thermal motions. This is 
expressed in terms of the local stability parameter Q being 
smaller than a critical value of order unity. 



< 



(1) 



is the radial velocity dispersion, (or the gas sound speed 
if it is larger), and k is the epicyclic frequency. We adopt 
K = \/3fl, where fl is the angular circular velocity. The 
factor -\/3 stands for a number between 1 and 2, depending 
on the shape of ^Ijr). For a thick disc, Qc — 0.68 
(iGoldreich & Lvnden-Belllll965h . 

The disc of mass Md is expected to develop transient 
elongated sheared features and to fragment into a few bound 
massive clumps of a typical mass 



and a typical radius 
i?c 0.52SRd , 



(2) 



(3) 



where 5 is the disc fraction of the total mass within the disc 
radius i?d. 



S = 



Ma 



(4) 



The fraction of disc mass in the clumps is crudely estimated 
to be in the range a ~ 0.1 — 0.4. Note that a higher 5 
corresponds to more massive clumps. 

An unstable disc is expected to self-regulate itself at 
Q Qc with a velocity dispersion-to-rotation ratio 



(5) 



The turbulence can be largely maintained 
by the gravitational encounters of the clumps themselves, 
which operate on a timescale <cnc — 2.1Q^a~^td, where 
id = — 50Myr is the disc dynamical time. With 

Q ~ 0.67 and a ~ 0.2, this matches the timescale for 
turbulence decay, idis — lAQ'^t^- The estimated effect 
of clump encounters is a lower limit to the total effect of 
the gravitational interactions involving all the components of 
the perturbed disc, which maintain the marginally unstable 
configuration. A higher S thus corresponds to a higher (Tr/V 
and a thicker disc. 

The same gravitational encounters and dynamical 
friction make the giant clumps migrate to the centre on a 
timescale 



t 



mig 



2.1Q^S'hd 



(6) 



and grow a bulge. The associated evacuation timescale 
for the entire disc mass is iovac — 10 .ba^^Q^S^^td- 
Additional comparable contributions to mass inflow in 
the disc are associated with angular-momentum transport 
outward, partly induced by the clump migration and partly 
due to torques involving the transient features. We noted that 
a higher 6 corresponds to more massive clumps. Then, With 
Q ~ Qc, a higher 6 or (Jr/V, namely more massive clumps, 
correspond to a more rapid migration. 

At high redshift, the continuous gas supply can maintain 
the disc unstable with giant clumps for seve ral Gigayears. 
The average cosmologic al gas accretion rate dNeistein et al.l 
120061: iGenel et al.ll2008l) is expected to be 



Here S is the surface density of the disc at radius r, cr^ 



M = 6 MiV^ (1 + z)2-25 /o.i6 Mq yr-1 (7) 

where M12 is the halo virial mass in units of lO^^M©, and 
/0.16 is the baryon fraction in units of 0.16. Based on the the 
effective deep pen etration seen in cosmological simulations 
dDekel et al.ll2009l) . the cosmological streams feed baryons 
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into massive galaxies on a timescale t^cc — 44r^^td, with 
r ~ 1 at 2 = 2 and varying from 2.5 to 0.4 between 
z ~ 9 and 1. The smooth component of the incoming 
streams, including small clumps, replenishes the evacuating 
disc, while the massive clumps associated with the incoming 
streams merge to the bulge. If the galaxy streams in which 
the massive clump fraction is less than average, the system 
is predicted to settle into a near steady state with 5 ~ 0.3, 
i.e., the input by streams and the transport from disc to 
bulge maintain a constant bulge-to-disc ratio near unity. This 
corresponds to comparable contributions of disc, bulge and 
dark matter to Aitot(^d)- 

On the other hand, in galaxies where the incoming 
streams are more clumpy than average, the merging of 
external massive clumps (possibly associated with dark 
matter haloes) are expected to grow a dominant spheroid 
with 5 < 0.25, and the dense clumps are capable of stirring 
up the turbulence in the disc to levels that stabilize the disc. 
DSC thus hypothesized that the dependence of the instability 
on the degree of dumpiness in the streams introduces a 
robust bimodality in the galaxy properties starting at z ^ 3. 

If the properties of the streams feeding a given galaxy 
vary in time, and in particular if the degree of dumpiness in 
these streams evolves, the galaxy may go through transitions 
from an unstable disc-dominated confi guration to a stable 
bulge-dominated state and vice versa (iMartig et al.ii2009l) . 
The instability can be responsible for efficient star formation 
while the stability may result in a significant suppression of 
star formation. However, after z ~ 2, the recovery from 
a bulge-dominated system back to an unstable disc takes 
several Gigayear and may never materialize. 

A systematic change in the stream properties is expected 
after z ~ 1 — 2. Then, the cosmological accretion 
rate becomes slower, the smooth cold streams no longer 
penetrate very effectively through the sho ck-heated media 
in massive haloes of ~ lO^^Mp, or hig her (Keres et al.'2 005l: 
Dekel & Birnboimll2006l: ICattaneo et al. 2006: Ocvirk etaLl 
2008h . and the input becomes dominated by stars rather 
than gas. Once the accretion cannot replenish the gas 
discs on a time scale comparable to the timescale for the 
discs to turn into clumps and stars, the galaxies gradually 
become dominated by old stars and eventually stable against 
axi-symmetric modes. This is because, unlike the dissipating 
turbulent gas, the stars maintain the velocity dispersion that 
they have acquired during their history. While the young 
stars are expected to have a velocity dispersion comparable 
to that of the gas and therefore to follow the gas into the 
bound clumps, the older stars are likely to have a higher 
(J. They therefore tend to form transient perturbations rather 
than accumulate in bound clumps, and eventually they join 
the stable component that just adds to the potential well and 
thus helps stabilizing the whole system. 

The axi-symmetr ic instability of a two-co mpo nent disc 
has be en studied by I Jog & SolomonI (1 19841) and iRafikovl 
(1200 lb . Denoting the velocity dispersions of stars and gas CTs 
and (Tg respectively (with the latter standing for the speed 
of sound if thermal pressure dominates), and defining Qg 
and Qg following eq. ^ separately for each component, 
the effective Q relevant for the instability of the combined 
system is approximately 



Table 1. Simulations details 



Comoving box size 
Number of DM particles 
DM mass resolution 
Max. resolution (proper) 
Min. mass of a stellar particle 



28.57 Mpc 
7 xlO" 
5.5 X W^Mq 
35-70 pc 



= 2( 



1 



is?' 



(8) 



where dgs = (Tg/ds and q is the dimensionless wave number 
q = fccTg/ K. The first term has to be slightly modified to take 
into account the dissipationless nature of t he stars but thi s 
correction makes only a small difference (IRafikovl 1200 Ih . 
The system is unstable once Q < 1. The most unstable 
wavelength corresponds to the q that minimizes Q; it lies 
between q = 1 for o-gs = 1 and q ~ a~^^ for Cgs <C 1. 
Note that with as > Ug, the stellar disc by itself may tend 
to be less unstable than the gas disc by itself, Qs > Qg, 
but through its contribution to the self-gravity that drives 
the instability, the stellar disc can help the gas component 
de-stabilize the disc. The combined system can be unstable 
even when each of the components has a Q value above 
unity. If, for example. Eg = Eg and (jgs ~ 0.5, the error 
made in Qc by ignoring the higher velocity dispersion of the 
stars is about 30%. If, on the other hand, CTgs ^ 1, the gas 
disc has to be treated on its own while the stars become part 
of the stabilizing component. 

The typical low-redshift discs are thus very different 
from the clumpy high-redshift discs. In a bulge-dominated, 
cold disc, supernova and stellar feedback may have an 
additional effect on the stabilization of the gas disc. 
The common low-redshift discs are expected to form 
predominantly in haloes below the threshold mass of 
^ lO^^MfT), and not necessarily by narrow streams 



( Birnboim & Dekel 


20031 Binnevl2004; Keres et al.l 120051; 


iDekel & Birnboim 


20061 iBirnboim et al.l 20071). Being gas 



poor, the low-redshift discs rarely develop axi-symmetric 
instabilities with Q < 1. They can evolve secularly 
through non-axi-symmetric modes of instability with Q ^ 1, 
associated with quiescent star formation. The turbulent 
high-redshift discs do not transform into today's thin discs. 
They end up most naturally in today's thick discs and SO 
galaxies, or in ellipticals through mergers. 

The above predictions of the simplified theoretical 
analysis are to be tested below using cosmological 
simulations. In parallel, the instability analysis will help us 
interpret the simulated results. 



3 THE SIMULATIONS 

The simulations that serve as the basis for the current 
analysis follow the evolution of three typical massive 
galaxies. They were performed completely independently, 
not tailored in any way to match the theoretical expectations 
or the observed SFGs. We selected dark-matter haloes with 
virial masses — lO^^Mo at z = 1, which end up 
as (3 — 4) X IO^^Mq today, comparable to and somewhat 
more massive t han the Milky Way. U sing the EPS statistic 
of halo growth (iNeistein et al.ll2006l) . the mean virial mass 
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of the major progenitor of such haloes is expected to be 
5 X lO^^M© at z ~ 2.3. The average input rate of gas 
into such a halo, assuming a universal baryonic fraction 
of /b = 0.165, is il/b ~ 6.6A/iV^(l + z)2-25 Mq yr^i, 
or Afb — 44AfQyr~^ for the haloes in question. These 
galaxies are thus at the low-mass end of the massive 
star-forming galaxies in the current surveys, in which the 
typical halo mas s is above and the ty pical SFR is 

~ 70Moyr-i jForster Schreiber et alj [20091) . This slight 
mismatch occurred because these simulations were not 
originally intended for the particular study reported in the 
current paper. 

The procedure for setting up the initial conditions of the 
simulations is as follows. We started with a low-resolution 
cosmological iV-body simulation in a comoving box of side 
20/i~^Mpc, using the following co smological parameter s 
suggested by the WMAP5 results (iKomatsu et alj |2009|) : 
fim = 0.27, nA = 0.73, 17b = 0.045, h = 0.7 and erg = 
0.82. At z = 1, we selected three haloes of — lO^^M© 
each, termed A, B and C. The selection was arbitrary, except 
that we deliberately avoided haloes that were involved in a 
major merger process at z ~ 1. This selection rejected less 
than 20% of the haloes of that mass and it had no explicit 
dependence on the merger history prior to z=l. One halo 
(C) happened to be in a relatively dense environment of 
neighboring haloes, and the other two (A and B) were in 
a "field" environment, namely part of a typical low-density 
filament of the cosmic web. These three haloes can therefore 
be regarded as a small fair sample of the haloes of that mass 
at that epoch. 

For each halo, we identified for re-simulation with high 
resolution the concentric sphere of a radius twice the virial 
radius, typically one comoving Mpc. This sphere was traced 
back to its Lagrangian volume at the in itial time, z = 60 , 
where we applied a zoom-in technique jKlvpin et al.ll200l1) 
to refine the fluctuations down to the desired resolution 
limit. Gas was added to the box following the dark matter 
distribution with the universal baryonic fraction /b = 0.165. 
We then re-simulated the whole box, with refined resolution 
only in the selected Lagrangian volume about the galaxy. 

Each simulation has a total of 7 x 10^ dark-matter 
particles, with four different masses. The high-resolution 
region is resolved with ^4x10^ dark-matter particles, each 
of mass 5.5 x IO^Mq. The particles representing stars have 
a minimum mass of 10^ Mq and their number increases as 
the simulation evolves; it reaches a value of 2 — 3 x 10^ 
stellar particles at z ~ 1.3. The integration of the gravity 
and the gas physics are performed on an adaptive mesh, 
where the maximum resolution is 35-70 pc in physical units 
at all times. This resolution is valid, for example, throughout 
the cold disc and dense clumps. The numbers characterizing 
each simulation are summarized in table[T] 

T he simulations were performed u sing the ART 
code jKravtsov et al.lll997l: [Kravtsovll2003h . following the 
evolution of the gravitating N-body system and the Eulerian 
gas dynamics using an Adaptive Refinement Tree. The 
code incorporates the standard physica l processes relevant 
for ga laxy formation, as described in ICeverino & KlvpinI 
(|2009|) . These include gas cooling by atomic Hydrogen and 
Helium, metals, and molecular Hydrogen, photoionization 
heating by a UV background with proper self-shielding, star 




Figure 1. Face-on gas surface density in galaxy C at 2 = 2.1. The image 
demonstrates violent disc fragmentation into transient features and bound 
clumps, resembling observed SFGs and theoretical expectations. The size of 
the image is 15 X 15 kpc. Surface density is in units of log(MQpc~^). For 
comparison, the surface density of molecular clouds in low-redshift galaxies 
is ~ 100 Mqpc-2, 

formation, metal enrichment, stellar mass loss, and stellar 
feedback. 

In our stochastic star-formation model, more than 90% 
of the stars form at temperatures well below lO^K and 
more than half the stars form below 300 K at gas densities 
higher than 10 cm^'^. New stellar particles are generated 
at every time step of the zero-level grid, which is about 
5 Myr at z = 2. We use a "constant" feedback model, 
in which the combined energy from stellar winds and 
supernova explosions is released as a constant heating rate 
over 40 Myr, the typical age of the lightest star that explodes 
as a type-II supernova. The heating rate due to feedback may 
or may not overcome the cooling rate, depending on the gas 
conditions in the star-form ing regions ( iDekel & Silklll986t 
ICeverino & KlvpinI |2009|) . We also include the effect of 
runaway stars by assigning a velocity kick of ^ 10 km to 
30% of the newly formed stellar particles. As a result, these 
stars can migrate ^ 100 pc away from the dense regions 
and explode as supernovae in regions with lower densities 
and longer cooling times. This enhances the efficiency of 
feedback in heating the surrounding gas. Finally, in order 
to mimic the self-shielding of galactic neutral hydrogen 
from the cosmological UV background, we assume that in 
regions where the gas density is higher than n = 0.1 cm~^, 
the ionizing flux shrinks to a negligible value (i.e., 5.9 x 
10^^ erg cm~^ Hz~^, the value of the pre-reionization 
UV background at z = 8). 



4 GLOBAL PROPERTIES AT HIGH REDSHIFT 

Our analysis here focuses on the three simulated galaxies at 
z ^ 2—2.5, when the Universe was roughly 3 Gyrold. These 
redshifts are in the middle of the redshift range 1.5 — 3 where 
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Table 2. Global properties of the host haloes at redshift 2 = 2.3. Radii 
are expressed in proper kpc, masses in units of 10^^ and velocities in 
km s~^. 



12 



Galaxy Rv Mv Vv 



A 
B 
C 



70 4 150 
68 3.5 140 
83 6.1 180 



the activity of disc instability and star formation is expected 
to be at its peak, both according to observations of SFGs and 
the predictions of DSC. Figure [T] presents the gas surface 
density in galaxy C at z = 2.1. The face-on view angle is 
defined by the angular momentum of the T < lO^K gas 
inside a sphere of 3 kpc. The image shows an extended disc 
of diameter ~ 9 kpc and the 70 pc resolution reveals that 
the disc is highly perturbed, with a few big, round clumps 
embedded in massive elongated tangential features. The 
elongated structures are transient features being disrupted by 
shear. The clumps are bound, with typical masses ~ IO^Mq 
and characteristic sizes ~ 1 kpc. About 10 — 20% of the 
total mass of the disc is in the clumps and sheared features. 
This galaxy, like the other two, also has a stellar bulge, of 
mass comparable to the overall disc mass. The very clumpy 
appearance of this high-redshift disc and its massive bulge, 
which repeats in the other two simulated galaxies, is revealed 
for the first time in cosmological simulations. A first visual 
inspection indicates that it resembles the clumpy appearance 
of the typical observed SFGs, and follows the theoretical 
expectations laid out by DSC, ^ 

Before we proceed to analyze the detailed properties 
of the simulated clumpy discs and bulges in comparison 
with the observed SFGs and the theoretical predictions, we 
address the main global properties of the dark-matter haloes 
and the stellar components in the simulated galaxies. We 
wish to verify in particular that they are consistent with 
the observed global properties of high-redshift star-forming 
galaxies, and their main scaling relations. 

4.1 Virial Properties of the Dark-Matter Haloes 

Table |2]displays the virial properties of the simulated haloes 
at z = 2.3. The virial radius and the corresponding mass 
Mv are defined to encompass a mean mass density of 180 
times the universal mean, and the virial velocity is = 
{GMy/R^y/^. Haloes A and B are similar, with ^ 4 x 
IO^^Mq and Vv — 150kms~^, while halo C is a bit more 
massive with AU = 6 x 10"Mq and K = 180 kms^^. 

Galaxy C has a growth history that differs from the other 
two. By z = 2.3, its halo has already assembled about half 
the mass that it will have at z = 1, while haloes A and B have 
assembled only a quarter of their z = 1 mass. The growth of 
C is thus faster prior to z = 2.3, and the growth of the other 
two is more rapid after z = 2.3. This reflects the fact that 
halo C was forming in an environment of higher density. The 
cosmological overdensity of matter in a concentric sphere 
of comoving radius about haloes A and B at 

z = 2.3 is only Sp/p ~ 0.07 and 0.02 respectively, while 
it is as high as Sp/p ~ 0.3 for halo C. This difference, 
representing the "cosmic variance", was set by our selection 
of haloes. We will see that the cosmic environment plays an 
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Figure 2. Scaling relations at z ~ 2.3 for our three simulated 
galaxies (blue squares) in comparison to observed galaxies (red circles). 
Observations are from the SINS survey (Forster Schreiber et al. 2003; 
ICresci et al.l l2009l). The top panel shows the Tully-Fisher relation, stellar 
mass versus maximum rotation velocity. The bottom panel shows the "main 
sequence" of star-forming galaxies, SFR versus stellar mass. The simulated 
galaxies seem to obey the observed scaling relations within the given scatter. 
They lie near the best-fit line of the Tully-Fisher relation, but they have SFR 
values about a factor of two lower than the average value for their stellar 
mass. 



important role in the formation of galaxies at the centres of 
these haloes, especially because the earlier formation time 
of halo C makes it denser in its central regions, both in dark 
matter and in gas. 



4.2 Scaling Relations at High Redshift 

Figure |2] refers to the two main scaling relations of SFGs, 
correlating the total stellar mass to SFR and to disc rotation 



High-redshift Clumpy Discs in Simulations 7 





Figure 3. A zoom-out gas surface-density maps showing tlie streams feeding galaxy A. Left: cold gas (T < 5 X 10* K). Right: hot gas (T > 3 X 10^ K). 
The box size is 160 X 160 kpc, covering the whole virial sphere. The color refers to log gas suiface density in units of H atoms cm~^ . Two major narrow 
streams carry the gas from well outside the virial radius to the inner ~ 20 kpc halo core, where they break into a multi-stream turbulent core before joining 
the inner disc of radius ~ 6 kpc, seen nearly edge-on at the box centre (mostly in white). 



velocity. The observational data is ext racted from the SINS 
survey dForster Schreiber et a n I2009h — a spectroscopic 
imaging survey in the near-infrared with SINFON WLT that 
focuses on SFGs at z 2. 

The top panel of figure |2]puts our simulated galaxies in 
the context of the TuUy-Fisher relation. The observational 
data is taken form the SINS sample of disc galaxies 
dCresci et al ] I2009h . Since the rotation curves in the 
simulated galaxies are rather flat outside the bulge, we adopt 
the average rotation velocity in the gaseous disc as our 
estimate for K„ax (see ^5.4| i. The simulated galaxies seem 
to follow quite well the observed Tully-Fisher relation. 

The bottom panel addresses the "main sequence" of star 
formation, namely the SFR as a function of total stellar mass. 
Given the scatter, the three simulated galaxies are consistent 
with the observed relation. The fact that they lie a factor 
of ^ 2 below the mean SFR at a given mass may be a 
random coincidence or it may reflect a certain deficiency 
in the simulated star-formation history. The stellar masses 
in A and B are at the level of 5% of their halo masses of ~ 
4 X 1O^^M0, which is quite typical for such haloes. The SFR 
in these galaxies are comparable to and only slightly smaller 
than the expected average gas accretion rate into haloes of 
similar masses, ^ 30 Mq yr^^. In galaxy C, with a halo only 
slightly more massive but a significantly denser environment 
and an earlier growth, the stellar mass is higher than average, 
and the SFR is higher than the average accretion rate. 

The agreement with the Tully-Fisher relation may 
indicate that the small systematic deviation in the SFR- A/, 
relation reflects a deficiency in the current SFR at z = 
2.3 rather than in stellar mass. It is thus possible that the 
simulated galaxies overproduce stars prior to z ~ 2.3, and 
then underestimate the gas fraction and the SFR at z ~ 2.3, 
while they reproduce the proper stellar mass. This may 
result, for example, from slightly underestimating the effect 



of supernova feedback at high redshifts. Indeed, the SFR, 
and how it is affected by stellar and supernova feedback, are 
naturally the most uncertain elements in our simulations, as 
these processes are implemented via crude recipes of subgrid 
physics. Despite this potential small deficiency, we conclude 
that the simulated galaxies do obey in a satisfactory way the 
main observed scaling relations, so their global properties 
can be considered as representative of typical SFGs of a 
similar mass at a similar time. 



5 CLUMPY DISCS AND BULGES IN STEADY STATE 
5.1 Feeding by Cold Streams 

Figure[3]shows a virial-scale view of the gas surface-density 
in halo A. The cold gas, T < 5 x lO'^K, flows along 
narrow streams that come from outside the virial radius 
and penetrate to the centre of the dark-matter halo. There, 
they blend in the inner ~ 20 kpc into a turbulent region 
and eventually settle in a clumpy disc of radius ~ 6 kpc. 
The shock-heated gas of T > 3 x lO'^K fifls much of the 
virial sphere. In the disc vicinity, the hot gas is also clumpy, 
perhaps indicating the presence of shocks at the interfaces 
between the streams and the central galaxy. 



5.2 Images of Gas, Stars and Dark Matter 

Figures 2116] show face-on and edge-on views of the galaxies 
in the inner ~'10% radius of the three haloes at z = 2.3. The 
maps display on the same scale the distribution of gas, SFR, 
stars and dark matter. 

Figure |4] refers to galaxy A. The gas surface density 
maps show a thick disc of radius ^ 6 kpc and thickness ~ 
1 kpc. The disc is highly fragmented showing several giant 
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Figure 4. Surface density maps of galaxy Aal z = 2.3. The side of each panel is 10 kpc. Shown are face-on views (left) and edge-on views (right). From 
top to bottom: gas, SFR, stars, and dark matter The units are log(M0pc~^) for the surface density and log( A/q yr~^kpc~'^) for the SFR. We see a disc of 
gas and young stars, a significant bulge extending to a very thick disc of old stars, and a rather spherical dark-matter halo. The disc is fragmented to elongated 
features and round dense clumps. Stars form at a high rate in the clumps and at the centre of the bulge. Most of the clumps are in-situ to the disc and are not 
associated with dark-matter concentrations. Two of the clumps are off-disc satellites. The color code was chosen to emphasize the clumps at the expense of 
saturating the highest densities. 
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Figure 5. Same as in figure|4]but for Galaxy B. This case shows the smallest disc of the sample. Beyond the high SFR in clumps, there is a relatively high SFR 
in other regions spread across the disc, some associated with the elongated transient features. The overall stellar population is rather smooth and concentrated 
in the bulge. Two clumps are external satellites with stars and dark matter. 
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Figure 6. Same as in figure|4]but for Galaxy C. This is the most massive and early forming galaxy in our sample. It has the highest surface densities of gas 
and stars and is therefore highly unstable. Stai' formation occurs in an extended clumpy ling. Two clumps are external, with stars and dark matter. 
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clumps and elongated transient clouds as in Figure [T] The 
SFR, as deduced from the stars younger than 10 Myr, reveals 
that a large fraction of the stars form in the disc clumps and 
in the centre of the bulge. The stars show a thick disc with 
embedded clumps, all associated with the gas, star-forming 
clumps, but the stellar distribution is dominated by a massive 
bulge, of a mass comparable to the disc mass. The bulge 
extends to a very thick disc, perhaps resembling a lenticular 
galaxy. The dark-matter distribution is rather spherical and 
smooth, with only two obvious dark-matter satellites, both 
off the disc plane. The disc clumps are thus all naked, 
not associated with small dark-matter haloes of their own. 
This is evidence for the in-situ nature of the disc clumps, 
indicating that they were originated in the disc itself rather 
than merged as small galaxies coming from the outside along 
the streams. The two off-plane satellite galaxies have haloes 
of masses 10^ — 10^ Mq and eccentric orbits that lead them 
to future mergers with the bulge while they never become 
part of the rotating disc. 

Figure |5] shows galaxy B. This galaxy, in the smallest 
halo of the three, has a disc radius of only ^ 3.5 kpc. The 
gaseous disc is strongly perturbed, with clumps embedded 
in elongated massive features that resemble spiral arms. The 
gas clumps here are less compact than in galaxy A, with 
a lower surface density extending to larger radii. The low 
stellar fraction in these clumps, compared to the clumps of 
galaxy A, reflects the transient nature of these features. The 
overall SFR of ~ 30 Mq yr^^ is similar to galaxy A but it 
is spread over the clumps and the broad elongated features 
covering a large fraction of the disc. The dark-matter map 
shows two satellites, coinciding with two clumps of gas and 
stars, of which one is actually in the disc plane. All the other 
clumps and elongated features are in-situ to the disc and not 
associated with dark-matter substructures. A non-negligible 
stellar bulge is seen, with a mass comparable to the disc 
mass, but it is less massive than in the other two cases. 

Figure |6] displays galaxy C, the most massive case that 
assembled most of its dark matter and accreted most of its 
baryons earlier than the other two. The gas disc is very 
clumpy and disturbed, despite the relatively low gas fraction 
of less than 10% in this galaxy. The instability is largely 
because the disc is more massive and more compact than 
in the other cases, with a high surface density in both gas 
and stars, ~ 2 x 10^ and 3 x IO'^Mqpc^^ respectively. 
The clumps and elongated features display a star-forming 
ring outside the bulge. This structure is remarkably similar 
to observ ations of star - formi ng outer rings in some z ^ 2 
galaxies dGenzel et alj l2008l) . Two off -plane satellites are 
seen in the stars, of which at least one is clearly associated 
with a dark matter satellite. The main features in the disc are 
not associated with dark-matter substructure. As in the other 
two cases, the bulge is as massive as the disc. 

In general, all three cases show a gravitationally 
unstable, extended, thick and irregular disc, with massive 
in-situ clumps and transient features, accompanied by 
a compact bulge of comparable mass, consistent with 
observations and the cosmological steady state predicted by 
theory. The mean gas surface density within the clumps 
is high, Eg 500 - 1000 Mgpc"^, compared to 

the typical values of ~ 100 Mqpc~^ in the "giant" 
molecular clouds in lo w-redshift galaxies jBolatto et al.l 
l2008t iHever et al.1 12008|) . This allows high, starburst-like. 



SFR surface densities of ~ 10 Mq yr ^ kpc as expecte d 
from the Schmidt-Kennicutt relation jKennicuttl Il998l) . 
Overall, 10 — 20% of the gas in the discs is in the 
form of clumps and sheared features. They are larger 
and much more massive than the star-forming molecular 
complexes in quiescent, low-redshift disc galaxies. This 
makes the gravitational interaction between the clumps 
and the transient features more intense, resulting in high 
turbulence and rapid mass flow to the center The surface 
density in the transient clouds is (100 — 5OO)M0pc~^ and 
they are embedded in a smoother inter-clump medium of 
surface density ~30 Mqpc~^. 

5.3 A Preliminary Comparison with Observed Images 

k 

In Fig. |2] we show side by side images of simulated 
and observed galaxies, as a sneak preview for a quantitative 
comparison to be provided elsewhere (Ceverino, Bournaud 
& Dekel 2009, in prep.). Edge-on and face-on images of 
simulated galaxy A at z = 2.3 are shown next to a 
"chain" and a "clump-cluster" galaxy respectively. Since the 
young stars contribute a significant fraction of the rest-frame 
U-band luminosity, the images of simulated stars younger 
than 100 Myr are compared to observed UV rest-frame 
images from the Hubble Ultra Deep Field (HUDF). 
The mean rest-frame UV (observed i band at z=2.3) 
mass-to-light ratio of the simulated stellar population is 
0. 1 Mq/Lq, ac cording to GALAXEV population synthesis 
code dBruzual & Chariot 2003) The pixel size and seeing are 
assumed to be 0.11 and 0.5 kpc, respectively. Galaxy A is 
shown at z = 2.3 because this is near where the clumpy 
appearance peaks in this particular simulation. The chosen 
UDF galaxies are at z = 3.3 and z = 1.3 respectively 
and they bracket our fiducial model at 2 = 2.3. This 
is fine for a very qualitative visual comparison given the 
observational evidence for a similar clumpy morphology 
of observed st ar-forming galaxies acr oss the redshift range 
z = 1 — 4 (e.g..lElme green et alj2007l ). The weak sensitivity 
to redshift is also consistent with the analytic prediction 
by DSC of a cosmological steady state. Therefore, for this 
very qualitative visual comparison, we did not make an 
effort to fully match redshifts, luminosities, masses and other 
galaxy properties, except for the overall similarity of disk 
size and the clumpy nature of the young stellar population. A 
visual comparison indicates that the general morphological 
features of the high-redshift clumpy SFGs are successfully 
reproduced by our simulations. At least in the specific cases 
shown, the similarity includes the number and sizes of 
clumps, their distribution throughout the system, and the 
contrast between the clumps and the inter-clump medium. 
This preliminary comparison is for illustrative purposes only 
— it is neither meant to provide a detailed modeling of the 
two specific galaxies, nor to present a statistically significant 
comparison of fair samples. 

The dominance of a few giant clumps is a common 
feature to the simulated and observed galaxies. For example, 
in galaxy A, about 40% of the young stars are in the in-situ 
disk clumps, and 10% in the bulge. The inter-clump medium, 
which is naturally more evident in the edge-on views, is 
the site for the other half of the SFR, which is associated 
with diffuse inter-clump UV light. This is consistent with 
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Figure 7. Typical images of observed and simulated clumpy disks at high redshifts. Left: UV rest-frame surface brightness of two galaxies from the Hubble 
Ultra Deep Field, representing the different appearances of the clumpy discs depending on the viewing angle: a "chain" galaxy (UDF 6478, z-band at 2 = 3.3) 
and a face-on "clump cluster" (UDF 9759, V-band at 2 = 1.2). The 3(t sensitivity in the images is 25.5 mag arsec^^ and the suiface biightness of the brightest 
clump (in black on the images) is around 21 mag arcsec"^. Right: The simulated galaxy A viewed at 2 = 2.3 with about the same resolution, using as a 
proxy for the UV emission the stars younger than 100 Myr The maximum and minimum surface densities are and 4 in log(MQpc~^). The color code is 
logarithmic and the same for the observations and the simulations. The pixel size and seeing are assumed to be 0.11 and 0.5 kpc, respectively. We do not 
aim here at modeling individual galaxies or making a valid statistical comparison between the simulated and observed moiphology. We pick the observed 
galaxies at redsh ifts that bracket our fid ucial simulation at 2 = 2.3, and because there is observational evidence for a similar morphology in the redshift range 
2 = 1 — 4 (e.g.. lElmegreen et alj200"i) . A visual comparison indicates that our simulations may reproduce the key morphological properties of high-redshift 
clumpy galaxies. 



the observation that ~^30% of the light in c l uster g alaxies 
is coming from clumps ( lElmegreen et alj l2005bl) . The 
clump sizes, their contrast with the rest of the disc, and 
their spatial distribution throughout the system, all seem 
strikingly similar between simulation and observation. In 
particular, some of the simulated galaxies show asymmetric 
configurations, resembling the observed "tadpole" galaxies. 
The mean stellar age within individual clumps, averaged 
over all the clumps of galaxy A, is ~'80 Myr The dispersion 
of stellar ages within a clump is comparable, ~70 Myr. This 
implies that the star-formation history within each clump is 
continuous over a period of 100 — 200 Myr. These stellar 
ages are similar to the age estimations from observations 
(lElmegreen et alj|2009l) . The mean stellar ages within the 
bulge is ~1 Gyr, with a spread of ^0.5 Gyr, in agreem ent 
with the observational estimates ( lElmegreen et al.ll2009l) . 

All three simulated galaxies host central stellar bulges 
with masses comparable to the disc mass, in agreement with 
the estim ated bulge-to-total mass ratio from the SFG s at high 
redshift dGenzel et al.1 120081: lElmegreen et al.ll2009l) . These 
stellar spheroids are rather compact, with radii (1 — 2) kpc. 



compatible with t he spheroids detected a t high redshifts 
(iKriek et al.ll2006l: Ivan Dokkum et aLll2008l: iBezanson et all 
120091) . 



5.4 Galaxy Properties vs Theory and Observation 

Table[3]lists global properties for the main components of the 
three galaxies: cold gas (T< 10^ K), stars and dark matter 
inside the disc radius, R^. This radius is estimated from a 
visual inspection of figures |4]|6] and the associated density 
profile of the stellar and gas disc, which are both truncated 
at a similar radius. The bulge radius i?b is read quite 
straightforwardly from the stellar density profile, which 
shows a well defined transition from a steep inner profile 
characteristic of a bulge and a shallower outer exponential 
profile with an exponential radius of 1 — 2 kpc, comparable 
to the bulge radius. The density profile of the gas disc is 
much shallower, with an exponential radius comparable to 
the outer disc radius R^- 

We use the kinematics of the star particles to distinguish 
the stellar disc from the spheroidal component. We assign 
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Table 3. Galaxy properties at 2 = 2.3. Radius in kpc and mass in IO^^Mq. 



Galaxy 


^disc 


^bulge 


Mgas 


Mstars 


Mdm 


Mbulgc 


Mdisc 


Mdisc /Mtotal 


Mgas /Mdisc 


Mbulgc/Mdisc 


■^clumps /-^disc 


A 


5.8 


2.0 


0.37 


1.73 


2.1 


1.0 


1.1 


0.26 


0.35 


0.93 


0.14 


B 


3.4 


1.0 


0.20 


1.27 


0.92 


0.81 


0.66 


0.27 


0.30 


1.22 


0.24 


C 


4.6 


1.5 


0.46 


11.44 


3.4 


6.6 


5.3 


0.35 


0.09 


1.25 


0.15 



Table 4. Star formation and baryonic fraction. SFR in Mq yr 
star-forming gas mass in 10^'' M0, and star-formation time in Myr 



Table 5. Kinematic disc properties. Velocities are in km s ^ . 



Galaxy SFR 


MsF_gas 


tSF 


Mgalaxy/Mv 


Galaxy Vgas 




(o'r/V)gas 


Vstars 


f^r, stars 


(o"r/V)stars 


A 27 


0.21 


78 


0.052 


A 180 


20 


0.11 


160 


33 


0.22 


B 25 


0.13 


52 


0.042 


B 180 


25 


0.14 


157 


34 


0.26 


C 70 


0.33 


47 


0.19 


C 380 


60 


0.16 


356 


107 


0.30 



a star particle to the disc only if the z-component of its 
angular momentum (parallel to the total galaxy angular 
momentum) is higher than a fraction /j of the maximum 
angular momentum for the same orbital energy, jmax = 
r, where \v\ is the magnitude of the particle velocity 
and r is its distance from the center. The distribution of 
jz/jmux among the simulated stars is bimodal, reflecting the 
division to distinct kinematic components of disc and bulge. 
We adopt as default fj — 0.7, which roughly marks the 
transition between the two components of this distribution. 
We thus define the bulge mass, Mbuigc, as the mass in stars 
inside i?d with jz/jmux > fj, and the disc mass, i\/diso 
as the mass of cold gas, Mgas, plus the mass in stars with 

jz/j max >/j. 

We find in the three galaxies values of S ^ 0.3 for 
the disc-to-total mass ratio (eq. (|4|i). These are compatible 
with the values predicted by DSC for unstable discs in 
cosmological steady state. The simulations thus reveal the 
predicted even distribution of mass between the disc, bulge 
and dark-matter components within the disc radius. 

The gas fractions in the discs of A and B are 35% 
and 30% respectively. This is significantly higher than in 
today's typical spirals, where the gas fraction is typically 
~ 10% or less. It is comparable, and perhaps on the 
low side of the crude observational estimates for typical 
high-redshift SFGs, both the indirect estimate s based on the 
Kenn icutt-Schmidt relation ( lErb et alj 120061: iGenzel et aTl 
I2OO6I) and the direct estimates based on CO observations 
dPaddi et all 12 004'; 'Greve et al.l 120051: iTacconi et all l2006l 
I2OO8I: iDaddi e t al. 2008, 2009). The gas fraction in galaxy 
C is lower than 10%, making it atypical. This may reflect an 
overproduction of stars too early in this simulation (see 0. 
The baryonic fractions within the disk radius in the three 
simulated galaxies, 0.5, 0,87, and 0.77, ar e consisten t with 
the observational estimates dForster Schr eiber et a n 120091 
Fig. 20). The fact that the gas fractions are lower than 
expected should be taken into account in the instability 
analysis, but it is encouraging to note that the clumpy 
appearance of the discs, and even the values of 5 as derived 
from the combined disc of gas plus stars, are consistent with 
the observations and the theoretical predictions. 

Assigning gas to clumps using a density threshold of 
n = 25 cm~'^ (or any value in the range (10 — 40) cm~'^), we 
obtain that the fraction of the disc in clumps is 0.15 — 0.25. 



Indeed, these fractions in clumps are consistent with the 
value of a ^ 0.2 adopted in the analysis of DSC. The mean 
clump mass in all cases is ^ IO^Mq, in agreement with 
the prediction, eq. (|2|i. The clump radius ranges from 0.2 
to 1 kpc, consistent with the pre-collapse radius predicted 
in eq. (O for galaxies of such mass, Rc ~ 0.8 kpc. The 
bound clumps, as opposed to the transient features, have 
rather sharp boundaries, and they are approximately in virial 
equilibrium. 

Table |4] lists galaxy properties that are related to star 
formation (discussed in section 14. 2t and baryonic fraction. 
The average SFR of 45 A/© yr^^ is in the ballpark of the 
expected mean gas accretion rates, eq. (|7]i, with the SFR 
in galaxy C slightly abov e the averag e . This is consistent 
with the main finding of iDekel et al.1 (l2009l) . based on a 
large hydrodynamical cosmological simulation, that a few 
cold streams carry almost all the accreted gas into the inner 
galaxy where it turns into stars in a rate comparable to the 
accretion rate. In order to estimate the timescale for star 
formation, or for gas consumption, we compute the mass of 
"star-forming gas", Mgp-gas, inside the disk radius R^. We 
define this as the gas with a density higher than 1 cm^^* and 
a temperature lower than lO"' K. It represents 60-70% of all 
the gas in the disc. The ratio between the star-forming gas 
and the SFR defines a SFR timescale of tsF = 50-80 Myr, 
which is somewhat sho rter than estimated from obs erved 
galaxies at high redshift ( iForster Schreiber et al.ll200"9l) . 

Two of the simulated galaxies show a low mass fraction 
of galactic baryons inside the dark-matter halo, which is 
25-32% of the cosmological baryonic fraction /b, probably 
reflecting gas outflows from the small building blocks of 
these galaxies. These values are in accord with the baryonic 
fractions that are estimated for massive, present-day galaxies 
(F. Prada, private communication). The baryonic fraction in 
galaxy C is high, and it may reflect an exceptionally high 
baryonic accretion-rate history. 

Table |5] summarizes the global kinematics properties 
of the gaseous and stellar disc: the rotational velocity, V, 
and the radial velocity dispersion, Cr- These quantities are 
the density-weighted mean values for the disc outside the 
bulge, r > i?b, where the rotation curve is roughly flat. The 
stellar disc is defined using the kinematic method described 
above. Galaxies A and B show similar overall kinematics. 
Their rotational velocity of F ~ 160 kms^^ is comparable 
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to that of today's discs in similar haloes. However, their 
velocity dispersion of cr^ — SOkms^^ is higher than in 
today's discs, reflecting a thick and irregular disc in which 
the pressure is dominated by turbulence. Within the clumps, 
the internal gas velocity dispersion is ~ 20 — 30 kms~^, 
and here too it is the dominant source of pressure. Galaxy C 
is again different; its rotational velocity and dispersion are 
significantly higher than in the other two cases. The <Jr/V 
ratio for the gas in the three cases is between 0.1 and 0.2, 
in agreement with the theoretical expectation, eq. Q, which 
is ar/V = 0.12 - 0.17 for 5 = 0.3 and Qc = 0.67 - 1.0. 
These values are compatible with the observed high-red shift 
discs (iGenzel et al .' 2006; F orster Schreiber et al.ll2009l) and 
are significantly higher than in low-redshift discs. 



5.5 Gravitational Instability 

Figure [8] displays face-on maps of the local effective 
Q parameter, computed following eq. ^ assuming a 
two-component disc of gas and stars. At each point in the 
disc, we used eq. ([T]l separately for each component to 
compute Qg and Qs from the corresponding values of S 
and cr,. and the common value of Vl. One can see that in 
all cases, a significant fraction of the disc is unstable with 
Q < 1. These fractions are 33%, 21% and 13% for A, B 
and C respectively. We also note that in galaxies A and B 
the most unstable regions define extended rings, similar to 
the rings of star fo rmation observed in some of the SFGs 
(iGenzel et al.ll2008h . 

The interpretation of the Q values obtained here from 
the simulated discs is not straightforward because they 
were measured form the highly perturbed discs while the 
linear analysis refers to the slightly perturbed disc at the 
onset of instability. We see that while the regions where 
Q < 1 naturally include all the dense clumps and transient 
perturbations, they also extend to regions that seem less 
perturbed, probably on their way to develop nonlinear 
perturbations. 

Galaxy A is the "most unstable" case, reflecting its 
high gas fraction and its low rotational velocity. Galaxy 
C is the most stable case, with a low gas fraction and a 
high rotational velocity. Still, galaxy C develops the most 
pronounced sheared features of the sample. About 60-70% 
of each disc have a Toomre parameter Q < 2, capable of 
developing high-order instabilities. It is interesting to note 
that these discs are highly unstable despite the dominance of 
the stellar component. As expected from eq. ([8]), the stellar 
component of the disc contributes to the local self-gravity 
that drives the instability, and helps destabilizing the disc. 



5.6 Evolution of Clumps and Spheroid Growth 

Figure|9]follows the formation of an individual clump in the 
disc of galaxy A near z ~ 2. It starts from a small-amplitude 
perturbation in a rather smooth region of the outer disc. 
During the initial 10 Myr, it is tangentially elongated due 
to shear. By 30 Myr it has turned into a high-contrast, 
round, bound and virialized clump. This is comparable to 
the dynamical time of this disc, = Rd/V ~ 31 Myr, or 
the free-fall time for a medium of n ~ 4 cm^-^. 

Some of the features seen in these images are transients. 




Figure 8. Instability analysis. Shown are maps of the effective Q parameter 
in the discs of galaxies A (top), B (middle) and C (bottom) overplotted 
over the gas surface density maps of figures |4]|1] and |6] The value of Q is 
computed from eq. JsJ assuming a two-component disc. The inner, solid 
contour refers to Q = 1, and the outer, dashed contour to Q = 2. In all 
cases, a significant fraction of the disc has Q < 1, so it is gravitationally 
unstable for axi-symmetric modes. The larger regions where 1 < Q < 2 
can be subject to higher modes of instability. 



where the strong shear prevents gravitational collapse. Such 
a transient can be seen in the bottom panels, rotating from 
left to top and starting to decay by the last snapshot. Similar 
elongated structures in the same galaxy are seen in the 
top panel of figure |4] They typically grow and decay on a 
dynamical timescale. 

Figure [TO] follows the evolution of an individual gas 
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Figure 9. Formation and evolution of a clump. Shown is a time sequence of disc A near z = 2. Each snapshot shows the face-on view of the maximum 3D 
gas density along the line of sight in units of Hydrogen atoms cm~'^. The size of each panel is 8kpc. A black square of size 1 kpc marks the position of 
one specific clump as it orbits the galaxy. The clump oiiginates as a small-amplitude perturbation in an unstable part of the disc and collapses to virialization 
within a dynamical time. 




Figure 10. Clump Migration into a bulge. Shown is a time sequence of disc A near z = 2. Each snapshot shows the face-on view of the maximum 3D gas 
density along the line of sight in units of Hydrogen atoms cm~^. The size of each panel is 6 kpc. A black square of size 1 kpc marks the position of one 
specific clump. The clump forms, virializes, and spirals in. It coalesces with the central bulge in 250 Myr, after completing about four rounds. 
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Figure 11. Evolution of galaxy B in steady state over cosmological times. Shown are face-on snap-shots of gas surface density at z = 2.8, 2.3, 1.6, 1.3 (from 
left to right) in units of log(M0pc~^). The box side is 10 kpc and the color code is the same in the three panels. These sequence demonstrates that this galaxy 
maintains a gravitationally unstable disc for at least 2-3 Gyr as it is intensely fed by cosmological cold streams. This steady-state ends near z = 1.3, when the 
disc is stabilized by the presence of a dominant bulge — an example of "morphological quenching". 



clump in galaxy A for a period of 250 Myr, from a 
pre-collapse cloud to its final coalescence into the bulge. 
At i = 0, we see an elongated broad perturbation. At 
t = 40 Myr, after about one dynamical time, the clump 
seems collapsed and virialized. As the clump rotates about 
the galaxy centre it spirals inward, while keeping its size 
and appearance roughly the same. By t = 250 Myr, after 
circling the centre about four times, the clump has coalesced 
with the central bulge. The migration time from the outer 
disc to the centre is ~ 8td, very similar to the prediction 
in eq. based on the gravitational interactions within the 
disc discussed in Sj2] The SFR inside the clumps starts of 
at a level of several Mq yr^^ and it gradually declines as 
the gas is substantially depleted. While most of the clump 
mass turns into stars before the migration is over, the clump 
is still gas rich when it coalesces with the bulge, leading to a 
compact bulge as in a wet merger. 

The migration of clumps to the centre contributes to the 
growth of the central spheroid. 10-20% of the disc mass 
flows into the bulge in ten dynamical times, in agreement 
with the estimated migration time by clump encounters, 
eq. (|6]l, and the fraction of the disc mass in the clumps, a ~ 
0.1 — 0.2. This is equivalent to one or two mergers of mass 
ratio 1:10 every 0.3 Gyr, which is more than expected on 
average for standard mergers of satellite clumps that com e 
with the streams from outside the galaxy jDekel et alj2009t) . 
Thus, the mass inflow in the disc contributes significantly 



to the early growth of spheroids. The bulge growth starts 
as soon as the disc becomes unstable, and it continues for 
several Gyr in a near steady state, as predicted by DSC. In 
this steady state, the continuous gas supply is replenishing 
the draining disc and the comparable bulge and disc keep 
growing in a similar pace. 



5.7 Steady State and Morphological Quenching 

Figure [TT| shows face-on snapshots of the gas disc of galaxy 
B from z = 2.8 to z = 1.3. During most of this 
period of 2.6 Gyr, the disc appears to be gravitationally 
unstable, with transient features and bound clumps, similar 
to its appearance at z = 2.3 as analyzed in detail 
above. This indicates a near steady state, as predicted 
by DSC. During this period, the gas supplied by the 
cosmological cold streams replenishes the gas that is 
being drained into the bulge through clump migration 
and angular-momentum transport. The value of (5 at z = 
2.8,2.3,1.6 is 0.23,0.27,0.25 respectively, demonstrating 
that this steady state is indeed characterized by a roughly 
even distribution of mass among the disc, bulge and dark 
matter halo components within the disc radius. A similar 
steady state is also seen in the other two galaxies. We 
conclude that the three galaxies simulated here represent the 
typical stream-fed galaxies observed as rapid star formers at 
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redshifts 3 to 1.5. A clumpy disc with a smaller bulge could 
appear within the first 0.5 Gyr after the onset of instability, 
perhaps following a period of stability due to a dense 
dark-matter halo, low accretion rate, or high turbulence. 

However, one can see in Fig.[TT]that by z = 1.3, galaxy 
B is no longer in the steady-state phase. The disc appears 
to be much smoother than before, showing only minor 
perturbations of much lower contrast. The disc fraction 
at this time is only 5 = 0.2, indicating that the disc 
is being stabilized by the presence of a dominant stellar 
bulge, constituting a significant fraction of the total stellar 
mass of Ah — lO^M©. As a result, the SFR in the 
disc is reduced to ~ 5M0yr~^, about a third of the 
average cosmological accretion rate at this epoch. This is 
an example of "morphological quen ching" of star format ion, 
as discussed and demonstrated in iMartig et alj ( l2009h . A 
similar phenomenon occurs in galaxy A. We note that 
the buildup of the stabilizing bulge in these galaxies is 
mostly due to mass inflow in the disc while in the case of 
iMartig et alj (|2009|) it is predominantly a result of external 
mergers. This demonstrates that morphological quenching 
can be obtained in different ways. 

We are not in a position to address the properties of 
the low-redshift descendents of our three simulated galaxies 
because the simulations were stopped at moderate redshifts, 
where they became too costly. The last snapshots available 
for galaxies A and B are at z = 1.3, and for galaxy C at 
z = 1.9. The appearance of a dominant stellar spheroid and 
a minor stable disc of stars and gas in galaxies B and A 
hints that these galaxies may end up as early-type galaxies 
with low SFR. On the other hand, as demonstrated by 
IMartig et al.l (l2009l) . such galaxies may eventually re-grow 
a sufficiently massive gas disc that becomes unstable again, 
and thus turn back into a blue, star-forming, late-type galaxy. 



6 ARTIFICIAL FRAGMENTATION 

If the Jeans length, Aj, is comparable to the resolution scale 
of a grid-based code, i.e., the smallest cell size A, one 
expects artificial fragmentation on that scale. This is because 
the code fails to properly follow the pressure gradients 
that are supposed to s upport a clump on that scale against 
gravitational collapse. iTruelove et alj ( Il997h have shown 
that the artificial fragmentation could be prevented if the 
Jeans length is always resolved by more than 4 resolution 
elements Q If we define a Jeans number by 



N, 



Aj 
A 



(9) 



it has to be kept above a critical value, = 4. The Truelove 
value was obtained based on simulations that lasted only a 
few free-fall times, so it is possible that a more stringent 
criterion may be required when the system is followed for 
longer periods, typical of cosmological simulations. 

Our current cosmological simulations do not resolve 
the thermal Jeans length of the coolest gas, T ^ 300K, 
so one must act to prevent artificial fragmentation. For this 

^ A similar criterion is necessary in SPH simulations, where the Jeans 
length has to be resol ved by more than a given number of SPH kernels 
iBate & BuAer3ll997l) . 





/ A 

Figure 12. Testing artificial fragmentation. Shown is the mass fraction 
in clumps (top panel), and the number of clumps (bottom panel) for 
simulations with different values of A^c- Convergence is indicated at 
Nc > 7, and artificial fragmentation becomes gradually more severe as 
Nc becomes smaller below Nc = 7. At the Truelove value, Nc = 4, the 
mass in clumps is larger than the base value at A^c 55 7 by a factor of 
~ 1.8 and the number of clumps is higher by a factor of ~ 2.3. The figure 
indicates that A^c 7 is the recommended value for a pressure floor that 
would safely eliminate artificial fragmentation. 



purpose we have implemented an artifici al pressure floor 
via a procedur e that is b ecoming stan dard (iMachacek et al.l 
2001; RobertsoiL&Kravtsov 2008; Ag ertz et al.l2009bl) . For 
gas with adiabatic index 7, the speed of sound is related to 
the pressure and density via Cg = jP/ p, so the Jeans length 
is given by 



A? 



(10) 
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In order to set the minimum iVj to a desired A^c, we thus set 
the pressure floor at 



P&oo 



(11) 



The pressure in the Euler equation is replaced by an effective 
pressure that is set to Pfloor when it would have otherwise 
obtained a lower value. We adopt 7 = 5/3. 

In order to verify the value of necessary for 
preventing artificial fragmentation in our cosmological 
discs, we re-ran the simulation of case B for the last 400 Myr 
prior to z = 2.3 with different values of the pressure floor, 
or Nc - The duration for this test is about 22 dynamical times, 
which should be sufficient for studying the fragmentation 
that develops on one or a few dynamical timescales. The 
corresponding values of Nc were chosen to be 2, 5, 7 and 
10, namely from half the Truelove value of A^c = 4 to more 
than twice above it. There is no point testing higher values 
of Nc because the corresponding Jeans length would have 
become comparable to the disk thickness and thus prevent 
the real disk fragmentation that we are after. 

Figures [12] shows the fraction of gas in clumps and the 
number of clumps as a function of A^c- The latter is measured 
as the number of density peaks with n > 25 cm^'^, and 
the trend is found to be insensitive to the exact threshold 
chosen. Below A^c = 7, we see a strong increase in both 
the mass in clumps and the clump number as a function 
of decreasing A^c- On the other hand, for A'c in the range 
7 — 10, we see a convergence of the two clump indicators 
into constant values. At the Truelove value, Nc ~ 4, both 
indicators have values about twice as large as their base 
values at A^c = 7 — 10. The similar behavior of the two 
indicators implies that the artificial fragmentation at small 
values of A^c is expressed in terms of the appearance of new 
clumps more than in the growth in mass of the real clumps 
that exist at large values of A'c. 

We conclude that in order to avoid 
artificial fragmentation, the Jeans length must be resolved by 
7 elements or more. We adopt A^c = 7 as our default case. 
With such a pressure floor, the clumps are likely to be real 
and their global properties are expected to be reliable. They 
are supported by turbulent motions in the disc and their size 
is set by the turbulence Jeans length. On the other hand, the 
detailed internal properties of the clumps, including density 
profile and substructure, are subject to the limited resolution. 
This internal structure can be resolved only in simulations of 
higher resolution. 



7 DISCUSSION AND SUMMARY 

We analyzed high-resolution AMR simulations of three 
galaxies in haloes of total mass ^ 5 x 1O^^M0 and 
baryonic mass ~ 5 x 10^°Mq at z ~ 2.3, selected quite 
randomly from a cosmological simulation. Each of these 
cases show a clumpy, thick and turbulent disc of gas and 
stars extending to (3 — 6) kpc in radius and accommodating 
a central stellar bulge of comparable mass within the inner 
(1 — 2) kpc. The discs show a very irregular morphology, 
with elongated transient features and a few bound massive 
clumps. The simulations reveal that these are naked clumps, 
not associated with dark matter haloes, and formed in-situ 



in the disc by gravitational instability. The simulated clumpy 
morphology is remarkably similar to the typical morphology 
of high-redshift SFGs as observed from different angles of 
view. 

The simulations i ndeed confirm many of the d etailed 
predictions outlined in lDekel. Sari & Ceverino! (l2009h based 
on a straightforward Toomre instability analysis in a 
cosmological context. They demonstrate that a large fraction 
of the massive galaxies at z ~ 2 — 3 are made of 
a gravitationally unstable disc with a stellar bulge of a 
comparable mass, corresponding to a disc-to-total mass ratio 
of (5 ~ 0.3 in eq. (HJi. A significant fraction of the disc has 
an effective value of Q ^ 1, typically defining an extended 
ring of maximum instability. The in-situ giant clumps are 
^ 1O®M0 each, namely on the order of one percent of the 
disc mass, in agreement with eq. (O. They involve 10-20% 
of the disc mass, as assumed in DSC. The clumps migrate 
to the disc centre in about 250 Myr, on the order of ten disc 
dynamical times, where they coalesce into a bulge. The mass 
inflow in the disc to the bulge contributes significantly to the 
spheroid growth, adding to the effect of mergers of external 
galaxies. 

In our simulations, the clumps survive bound until 
they merge into the bulge. It has been argued that 
in reality the clumps are disrupted on a dynamical 
timescale by radiative feedback from newly formed stars 
( iMurrav. Ouataert & Thompson .2009). In an associated 
paper (Krumholz & Dekel, in preparation), it is shown 
that the ejected fraction of gas from the giant clumps 
in massive disks at high redshift is limited to less than 
10% provided that the SFR efficiency is ^ 0.01, i.e. it 
obeys t he Schmidt-Kenn icutt relation between SFR and gas 
densitv pCennicuttll998h. which is also favored theore tically 
dKrumholz & McKed l2005t iKrumholz&T^ |2007|) . The 
clumps could be disrupted only if the SFR efficiency was 
larger than its local value by an order of magnitude. 

The SFR at z = 2.3 is at the level of ~ 40 Mq yr^\ 
which is in the ball park of the expected cosmological gas 
accretion rate along cold streams into haloes of comparable 
masses foekel et al. 2009.) . This corresponds to more than 
lOOM0yr~^ in the massive SFGs that reside in haloes 
more massive than IO^^Mq. We recall, though, that the 
simulated SFR at z ~ 2.3 may be slightly underestimated 
due to the limitations of the sub-grid recipes for star 
formation and feedback. Given these limitations, about half 
the star formation seems to occur in the disc clumps. The 
star formation within each clump extends to ^ 200 Myr, 
or several disc dynamical times, which is comparable to 
the timescale for clump migration. If the star-formation 
efficienc y is ~ 1 — 2% as in the molecular clouds of today's 
galaxies dKrumholz et alj2009l) . then it must occur in denser 
sub-clumps (DSC), which are not resolved in our current 
simulations. 

Our simulations confirm that the disc is continuously 
supplied by fresh material through cold streams that 
penetrate the hot halo along the dark-matter filaments of the 
cosmic web. In a companion paper, we use the simulations 
to predict the observational appe arance of these col d streams 
as extended Lyman-alpha blobs dGoerdt et alj2009l) . Despite 
the very high resolution, the streams in our three galaxies 
are rather smooth; the fraction of the incoming stream mass 
in clumps with densities higher than 0.1 cm^'^ ranges from 
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7% to 25%, i n compariso n with the average value of 30% 
estimated by iDekel et aTl (|2009|) for massive galaxies at 
similar redshifts. As argued by DSC, this smoothness of the 
incoming mass is a key for maintaining the disc sufficiently 
dense and cold without growing a stabilizing bulge, such 
that the disc is kept gravitationally unstable. The continuous 
supply of gas keeps this configuration of an unstable disc 
and a comparable bulge in a near steady state that lasts for 
2-3 Gyr. 

At least two of the galaxies simulated here become 
rather stable before z = 1; they become dominated by a 
stellar spheroid that stabilizes the remaining disc, prevents 
its fragmentation into giant clumps, and suppresses star 
formation to quiescent levels. This is a demonstration 
of the morpho l ogical- quenching mechanism discussed 



by Martigetal. ( 20091) and seen in the simulation of 



lAgertz etal.l(l2009bl)" 

Two of the simulated galaxies show a mass fraction of 
galactic baryons inside the dark-matter halo which is ~ 30% 
of the cosmological baryonic fraction, while the third galaxy 
has a baryon fraction comparable to the universal value. 
The comparison with the low baryonic fractions indicated 
from observations at low redshift (e.g., F. Prada, private 
communication) is beyond the scope of this paper, because 
our simulations are currently limited to high redshift. We 
anticipate that low baryonic fractions can be obtained if the 
baryon accretion is suppressed at low redshifts while the 
dark-matter haloes keep growing. 

While the current simulations are state of the art, they 
are still limited in many ways. One limiting issue is the 
relatively low gas fractions of 8 — 35% in the simulated 
discs compared to the values of 30 — 60% indicated for 
the observed galaxies. Although the observational estimates 
are indirect and therefore uncertain, it is possible that our 
current simulations slightly overestimate the star formation 
rate at earlier times, possibly as a result of underestimating 
the suppression of SFR by stellar and supernova feedback. 
This may lead to an underestimate of the gas fraction and 
SFR at z — 2.3. Fortunately, this possible limitation does 
not prevent the robust appearance of strong instability in the 
high-redshift discs at z ^ 1.5, so the simulations do produce 
galaxies that mimic the main features of the observed SFGs 
and the theoretical predictions. 

Another caveat is associated with the mandatory 
incorporation of a pressure floor. While it serves the 
purpose of preventing artificial fragmentation, it necessarily 
smoothes the thermal-pressure gradients in high density 
regions. Only simulations with a resolution better than 10 pc 
could help us directly evaluate the side effects that may be 
introduced by this pressure floor on small scales. However, 
the tests presented in fj6] demonstrate that the global 
properties of the clumps are stable against variations of the 
pressure floor in the regime where artificial fragmentation 
is suppressed. This is because the gravitational collapse of 
the clumps is balanced by velocity dispersion rather than by 
thermal pressure gradients. The velocity associated with the 
artificial pressure is ~ 8kms^^ compared to the turbulent 
velocity of 20 — 30kms^^. In practice, over most of the 
simulated unstable disks, the Jeans length associated with 
the overall pressure is significantly larger than the artificial 
minimum imposed at 7 cell sizes. 

Our main findings from the clumpy galaxies at 



high redshif t are co nsistent with the parallel findings of 
lAgertz et alj (l2009bh . who used a simulation with a similar 
resolution but a different AMR code. They also see clump 
formation in massive spiral arms driven by disc instability. 
However, they find about twice as many bound clumps, 
ranging to smaller masses and preferentially populating 
the outer disc, which seem not to resemble the observed 
SFGs as well as our simulations. Furthermore, they detect 
significant small-scale fragmentation in the initially smooth, 
cold streams, where we see no such effect. This is probably 
a result of some degree of artificial fragmentation in 
lAgertz et alj (l2009bh due to the lower pressure floor adopted 
in their simulation, (f|6ll. When resolving the Jeans length 
with only 4 resolution elements, following the Truelove 
criterion, we find that artificial fragmentation tends to occur 
in regions where the gas velocity dispersion is low such that 
the Jeans length is strongly affected by thermal pressure, 
as in the coherent cold streams. On the other hand, when 
we resolve the Jeans length by 7 resolution elements, we 
prevent such artificial fragmentation in our simulations. Still, 
the man y qualitative similarities between the findings of 
Agertz et al. (2009b) and ours indicate that we both capture 
the real key features of high-redshift SFGs. 

We conclude that zoom-in simulations of galaxies 
drawn quite randomly from a cosmological volume 
reproduce clumpy discs and bulges similar to the observed 
star-forming galaxies at high redshift. The simulations 
demonstrate the validity of the simple theoretical picture 
where the continuous gas supply by cold streams drives 
a self-regulated gravitational instability in a cosmological 
steady state that tends to last for several Gyr. 
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